Neutral pion photoproduction off protons in the energy range 0.3 GeV < E 7 < 3 GeV 
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Single 7T° photoproduction has been studied with the CB-ELSA experiment at Bonn using tagged 
photon energies between 0.3 and 3.0 GeV. The experimental setup covers a very large solid angle 
of ^98% of 4-7T. Differential cross sections da/dil have been measured. Complicated structures in 
the angular distributions indicate a variety of different resonances being produced in the s channel 
intermediate state 7p — > N*(A*) — > p7r°. A combined analysis including the data presented in this 
letter along with other data sets reveals contributions from known resonances and evidence for a 
new resonance N(2070)Dis. 

PACS numbers: 14.20 



The study of baryon resonances, of their masses, 
widths, and decay modes, is a cornerstone to improve 
our understanding of the internal structure of nucleons. 
Apart from the recently suggested exotic pentaquark 
state 9 + (1540) 1], the pattern of known baryon reso- 
nances below masses of 1.8 GeV listed by the Particle 
Data Group (PDG) |2| is generally well reproduced 
by present quark model calculations 0- d ■ Above 
this mass, models predict many more resonances than 
have been seen experimentally. Moreover, the mass 
predictions of models using one-gluon exchange [3j or 
instanton-induced interactions |5| no longer agree. Lat- 
tice gauge calculations reproduce the masses of ground- 
state baryons rather well jfjj and even first excited states 
have been simulated However the structure of 

baryons is still far from being understood. 

Results of elastic and inelastic 7rN scattering experi- 
ments provide the largest fraction of our knowledge on 
the excitation spectra of the nucleon and the A(1232)P33. 
Most excited states are also reached in photonuclear reac- 
tions, yielding information on their photocouplings. The 
investigation of the neutral pion channel is particularly 
interesting since background terms from direct produc- 



tion off the pion cloud or from production via i-channcl 
exchange of pions are suppressed. 

Experimental information on ir° photoproduction is 
sparse at high photon energies. Our existing knowledge 
is condensed into the parameterizations of SAID, a par- 
tial wave analysis program run at GWU @, and of 
MAID, the Mainz unitarity isobar model for pion photo- 
and electroproduction . Most data sets cover the low- 
energy region or date back to the 1970s and are limited 
in solid angle coverage and accuracy. Here, we just quote 
the recent JLab data on single 7r° photoproduction, cov- 
ering the photon energy range from 434 to 1742 MeV 
and an angular range which varies from —0.89 to 0.15 
in cos6* cm at the lowest, and from —0.79 to 0.89 at the 
highest energy [llj . 

In this letter we report on a new measurement of sin- 
gle 7T° photoproduction off protons covering the energy 
range from the A(1232)P33 mass region up to an invari- 
ant mass of 2.54 GeV/c 2 . The experiment (Fig. [TJ was 
carried out at the tagged photon beam of the ELectron 
Stretcher Accelerator (ELSA) at the University of Bonn. 
Electrons were extracted at energies of 1.4 and 3.2 GeV, 
covering photon energies from 0.3 to 3.0 GeV, with a typ- 
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ical intensity of 1-3 xlO 6 tagged photons/s. The photon 
beam hit a liquid H2 target of 5 cm length and 3 cm di- 
ameter. Charged reaction products were detected in a 
three-layer scintillating fiber (scifi) detector. One of the 
layers was straight, the fibers of the other two layers en- 
circled the target with ±25° with respect to the first 
layer. The intersection point of a proton could be re- 
constructed for polar angles from 15° to 165°. Charged 
particles and photons were detected in the Crystal Barrel 
detector |l2J. It consisted of 1380 CsI(Tl) crystals with 
photodiode readout covering 98% of An. A total absorp- 
tion photon detector (oil Cerenkov counter, not shown 
in Fig. 1) placed downstream was used to determine the 
photon flux. 

The coincidence between tagger and scifi detector pro- 
vided the first-level trigger of the experiment. From the 
hit pattern in the Crystal Barrel detector, a fast clus- 
ter logic determined the number of 'particles' defined by 
clusters of contiguous crystals with individual energy de- 
posits above 15MeV. A second level trigger was gener- 
ated for events with two or more 'particles' in the cluster 
logic. In the data analysis, clusters with two local max- 
ima were split into two 'particles' sharing the total energy 
deposit. The offline threshold for accepted particles was 
set to 20 MeV. The proton kinetic energy had to exceed 
35 MeV to traverse the inner two scifi layers and to pro- 
duce a trigger. A proton energy of 90 MeV was needed to 
reach the barrel calorimeter and to deposit the minimum 
cluster energy of 20 MeV. 

In the first step of data reduction, we required both 
photons to be detected in the Crystal Barrel and that not 
more than one candidate for the proton was found in bar- 
rel and/or inner detector. Thus, we retained events with 
two or three 'particles'. Proton candidates were identi- 
fied by the geometrical relation of impact points in the 
scifi and in the barrel. In the further analysis, protons 
were not explicitely used but treated as missing parti- 
cles. The remaining two 'particles' were treated as pho- 
tons. These events were subjected to kinematic fits im- 
posing energy and momentum conservation by fitting the 
p(7! 77)Pmissing hypothesis using one constraint (1 C) and 
the p(7, 7r°)p m i ss i ng hypothesis (2C) for E 1 <1.3GeV. 

Fig. 2 shows the 77 invariant mass spectrum after a 
10~ 4 confidence level cut in the 1 C kinematic fit. The 
7T° and 77 signals are observed above a small residual back- 
ground. There are 2.6 x 10 6 events due to 7p — > p7r°. The 
data is listed numerically in the Durham data base [l3j |. 
The photoproduction of rj mesons is discussed in a sub- 
sequent letter |l4j . 

The detector acceptance was determined from 
GEANT-based Monte Carlo simulations. It vanishes for 
forward protons leaving the Crystal Barrel through the 
forward hole, and for protons with very low lab momenta. 

The absolute normalization was derived from the mea- 
sured photon flux and by a comparison of our data with 
SAID. For each energy channel in the tagger, the an- 
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Figure 1: 
Experimental 
setup at ELSA in 
Bonn 

Figure 2: Spec- 
trum of invariant 
77 masses on a log- 
arithmic scale. 



gular distributions were divided into 20 bins in cos 9 cm 
where 9 cm denotes the emission angle of the meson in 
the center-of-mass system with respect to the initial pho- 
ton direction. For the low energy data (below 1.3 GeV) 
the resulting distributions were fitted to SAID by the x 2 
method for each energy channel. In general, the shapes 
of the distributions are in excellent agreement. The ab- 
solute flux normalization was determined from the SAID 
data base. For the high-energy data, this method could 
have been applied up to about 2 GeV. Above, experi- 
mental and SAID angular distributions deviate signifi- 
cantly. Below 2 GeV, a comparison of our angular dis- 
tributions and SAID defines the absolute normalization 
of our photon flux, with one scaling factor over the full 
energy range. The uncertainty of the normalization is 
estimated to 5% for photon energies below 1.3 GeV and 
to 15% above. 

The differential cross sections da/d£l are shown in 
Fig. |21 Statistical and systematic errors are added 
quadratically. The latter were calculated from Monte 
Carlo simulations of different experimental uncertainties, 
a shift of the photon beam by 3 mm with respect to the 
target center, a target shift of ±1.5 mm in beam direc- 
tion and a change in the thickness of material between 
barrel and inner detector corresponding to 1 mm of cap- 
ton foil. The relative reconstruction efficiency error was 
estimated to 5%. As the angular distributions are not 
influenced by the flux errors, these are not included in 
Fig- El The dark gray line shows the result of the partial 
wave analysis described below. The total cross section in 
Fig. 0] was obtained by integration of the measured cross 
section and extrapolation to the non-measured region by 
using the results of the partial wave analysis (P WA) . 

The data was subjected to a coupled-channel par- 
tial wave analysis within the framework of an isobar 
model. Resonances are described by Breit-Wigner am- 
plitudes, except for strongly overlapping states which are 
described in K-matrix formalism. Angular distributions 
are calculated in the spin-orbital decomposition formal- 
ism [HOH. The PWA method is described in [HI- 
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Figure 3: Differential cross sections for 7p — ■> p7r° (E-y bins in MeV). COS 8 cm 
The measured data points are shown as black squares, the PWA result as a dark, 
the SAID parameterization as light gray line. At large _B 7 , SAID reproduces the 
-1 -0.5 0.5 1 -0.5 0.5 1 very forward data points, yet the agreement over the full angular range is worse. 



Along with the CB-ELSA data on photoproduction 
of 7T° and r] mesons as discussed here and in the subse- 
quent letter Q, several other data sets served as in- 
put for the analysis: The TAPS low-energy data on 
7P — > pr] 18j to fix the 77 threshold; beam-asymmetry 
measurements of tt° [ljj, |2fJ and 77 [ljj, [21j photoproduc- 
tion were used to provide additional constraints from po- 
larization observables. The beam-asymmetry data cov- 
ers the mass range from 1200 to 2000 MeV. We include 
data on 7p — > n7r + [2^ covering the mass range from 1500 
to 2150 MeV thus exploiting the different decay couplings 
of N* and A* resonances to nir + and p7r°. 

The x 2 °f the final fit values for the different data sets 
are given m 14]. The masses, widths, helicity ratios, 
and contributions of different resonances to the cross sec- 
tion for 7T° photoproduction from the final solution are 
given in Tab. HJ We consider contributions to the to- 
tal cross section with fractional yields (column 5) above 
1-2% as established from this analysis. The sum of frac- 



tional yields is not equal to one due to interferences. 

In the total cross section clear peaks are observed for 
the first, second, and third resonance region. A fur- 
ther broad enhancement can be seen at about 1900 MeV. 
The first resonance region is of course mainly due 
to A(1232)P 33 ; the N(1440)P U Roper resonance pro- 
vides a small contribution of about 3-4% compared to 
A(1232)P33. The Ai/ 2 /A 3 / 2 helicity ratio of excitation of 
A(1232)P33 is determined to be 0.52 ± 0.03 which agrees 
favorably with the PDG average 0.53 ± 0.03. 

The second resonance region receives approximately 
equal contributions from N(1535)Sn and N(1520)Di 3 . 
The helicity ratio \A 1 / 2 /A 3 / 2 | of the excitation of the 
N(1520)D 13 is determined to be small (< 0.18). From the 
couplings given by the PDG we calculate —(0.14 ± 0.06); 
Mukhopadhyay used an effective Lagrangian approach to 
determine it to be -(0.4 ± 0.1 ± 0.1) 23]. Within ~ 2a 
there is consistency between these values. 

The third bump in the data is made up out of three ma- 
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Table I: Masses, widths and resonance couplings, this analy- 
sis. 
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* K matrix fit, pole position of the scattering amplitude in the 



complex plane, fraction for the total K-matrix contribution 
° This contribution includes non-resonant background. 
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Figure 4: Total cross section for 7p — » p7r° obtained by in- 
tegration of angular distributions and extrapolation into for- 
ward and backward regions using of PWA result. The dark 
gray line represents the result of the PWA, the light gray line 
represents SAID. 

jor contributions. The A(1700)D33 provides the largest 
fraction (~ 35%) of the peak, followed by N(1680)Fi 5 
(~ 25%) and N(1650)Sn (~ 20%) as extracted from 
K-matrix parameterization; A(1620)S3i (~ 7%) and 
N(1720)Pi3 (^6%) are also observed. In the fourth reso- 
nance region we identify A(1950)F37 contributing ^41% 
to the enhancement and A(1920)P33 with ^35%. 

The fit also suggests the presence of A(1905)F35 and 
A(1940)D33. Additionally, weak evidence is found for 
the new resonance N(2070)Dis reported in in the 
p7r° decay mode (~ 9%). Here we mention only that 
N(2070)Di5 is predicted by quark models as being one of 



the missing states. In |3j, a doublet of ND15 resonances 
is predicted at 2080 and 2095 MeV, respectively, where 
the lower (and only the lower) mass state has a large 
coupling to N?7. A similar mass pattern is found in |5j. 
The high-energy region is dominated by p-ui exchange in 
the t channel as can be seen by the forward peaking in 
the differential cross sections. 

We have reported new data on the photoproduction 
of neutral pions off protons over the full baryon reso- 
nance region. The data extends the existing data base in 
larger solid angle coverage and in energy. A partial wave 
analysis of the data presented here and elsewhere 0, 
decomposes the data into contributions from a large num- 
ber of N* and A* resonances. 
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